Statics and dynamics of a cylindrical droplet under an external body force 
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We study the rolling and sliding motion of droplets on a corrugated substrate by Molecular 
Dynamics simulations. Droplets are driven by an external body force (gravity) and we investigate 
the velocity profile and dissipation mechanisms in the steady state. The cylindrical geometry allows 
us to consider a large range of droplet sizes. The velocity of small droplets with a large contact 
angle is dominated by the friction at the substrate and the velocity of the center of mass scales 
like the square root of the droplet size. For large droplets or small contact angles, however, viscous 
dissipation of the flow inside the volume of the droplet dictates the center of mass velocity that 
scales linearly with the size. We derive a simple analytical description predicting the dependence 
of the center of mass velocity on droplet size and the slip length at the substrate. In the limit of 
vanishing droplet velocity we quantitatively compare our simulation results to the predictions and 
good agreement without adjustable parameters is found. 



I. INTRODUCTION 

Understanding and controlling the motion of droplets 
is crucial for many practical applications including ink jet 
printing or microfluidic devices^ The motion of droplets 
can be driven in various ways: For instance, if a tem- 
perature gradient is exerted on the fluid the droplet will 
move because of the Marangoni forces caused by sur- 
face tension gradients;^ Alternatively, a temperature^ or 
wettability gradient^i^ can be applied at the substrate. 
Forces can also be exerted onto a droplet by immersing 
it into a sheared, viscous fluidj2i 

Instead of transferring forces via the droplet 's sur- 
face one can apply a body force. Examples are falling 
drops on an inclined plane due to gravit y^i^i^° or the 
motion of droplets on a rotating disk due to centrifu- 
gal forces. In these cases much attention has been de- 
voted to the onset of motion and the changes of the 
shape of the droplets and their instabilities in response to 
strong body forces fi^'i^'i^iii^'i^ Surprisingly, the steady- 
state motion of droplets under these conditions has re- 
ceived little attention. ^S. However, this linear regime, 
where the droplet shape resembles its equilibrium form, 
is pertinent to small droplets as their occur in nanoflu- 
idic devices. In the case of small droplets one can ask the 
question by which mechanism the energy that is imparted 
onto the drop via the external force is dissipated. The 
dominating mechanism is expected to depend on the size 
of the droplet^^ and it dictates the steady-state velocity 
of the droplet. 

Similar dissipation mechanisms also act when a droplet 
spreads on a wettable substrat a^''''-^^i-^^'^°i^^'^^i^^i^'^ but 
the study of the steady-state properties is computation- 
ally much more convenient because one can average the 
properties (like the velocity profiles inside the droplet 
shown in Fig. [1]) along the trajectory instead of averag- 
ing over many realizations of the spreading process. 

Here we study the steady-state of cylindrical droplets 
under the influence of a weak body force and determine 
the flow profiles inside of the droplet (see Fig. [T|) as well 
as the relation between the velocity of the droplet 's cen- 



ter of mass and its size in the linear regime. We derive 
a simple analytical approximation for the velocity inside 
the droplet, which is quantitatively compared to the re- 
sults of extensive Molecular Dynamics simulations of a 
coarse-grain polymer model. 

The manuscript is organized as follows: In Sec. |TT] we 
describe our simulation model and technique. Then, in 
Sec. mil we compute properties of the bulk liquid such as 
surface tension and viscosity. Sees. IIVI and |V] comprise 
the analysis of static properties, such as contact angle and 
adhesion energy, and dynamics properties of droplets, re- 
spectively. The manuscript closes with a brief discussion. 
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FIG. 1: Velocity field inside a droplet of contact angle 6e ~ 
130°. The inset represents the viscous energy dissipation per 
unit volume of fluid as a function of the height. 
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II. THE MODEL 

Being interested in the universal properties of liq- 
uid droplets we utilize a coarse-grain model of a poly- 
mer liquidi^i2£i^iS^ In experiments as well as simulation 
studies, polymer liquids are particularly suitable because 
the vapor pressure is very low and evaporation effects 
can be disregarded, Moreover, the chain length allows 
to control the viscosity without changing the molecu- 
lar interactions. Our coarse-grained polymer model in- 
corporates only the three relevant interactions necessary 
to bring about the universal properties of dense poly- 
mer liquids: connectivity of the segments along the ex- 
tended macromolecule, excluded volume of the segments 
on short distances and an attractive interaction at inter- 
mediate distances to make the polymers condense into a 
liquid. Within the framework of such a coarse-grained 
model, an interaction center - a segment - is comprised 
of a small number of monomeric repeat units. 

Our model for the polymer droplets is similar to the 
one of Pastorino et al. for a polymer melt interacting 
with a polymer brush. The potential between neigh- 
boring monomers of a polymer is given by the Finite 
Extensible Nonhnear Elastic (FENE) potential: 
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with Rq = 1.5 (7 and fc = 30 e/cr^. The monomers addi- 
tionally interact with a 6-12 Lennard- Jones potential. 
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where the cutoff distance is chosen to be rc = 2 X 2i/6 a. 
The Lennard- Jones parameters are fixed to unity, e = 
1 and (7=1. In all the simulation the length of the 
polymers is fixed to Np — 10 monomers. 

The solid, corrugated substrate is modeled by two lay- 
ers of an fee lattice of density ps — 2 a^^ interacting with 
the monomers via a Lennard- Jones potential. The length 
scale of the potential is fixed to CTs — 0.75 cr, while the 
strength is varied between = 0.2 — 1.0 e. The cutoff 
distance is identical to the range of the interactions in 
the fluid. This attraction between the polymer segments 
and the substrate dictates both static properties (e.g., 
contact angle) and dynamic properties (e.g., slip length). 

We use a dissipative particle dynamics (DPD) ther- 
mostat in our model. The DPD thermostat has the ad- 
vantage of acting locally on the system while conserving 
the momentum. Hence the occurrence of non-physical 
currents in the fluid is avoidedi^SiSii^^ With a DPD ther- 
mostat two new force are added to the microcanonical 
equation of motion, 



(3) 



where m denotes a segment's mass and f its acceleration. 

and F^ are the dissipative and fluctuating forces act- 
ing on particle i due to the interaction with particle j, 
respectively. In order to conserve momentum these terms 
are applied to pairs of particles and they take the form: 



F^ = C^nin 



(4) 
(5) 



where l.y — Vij/vij is the unit vector pointing from par- 
ticle j to particle i. The damping coefficient, 7, and the 
strength of the noise, C, obey the fluctuation dissipation 
theorem, = 2kBT^. We fix 7 = 0.5 in all our simula- 
tions. The weight functions are defined as. 



(1 - t/tcY r < rc 
r > Tc 



(6) 



and Oij in Eq. ([5]) is a random variable of vanishing av- 
erage and second moment: 



= 



(7) 
(8) 



We use a uniform distribution for 9 because it was shown 
that uniform random numbers give the same results as 
Gaussian random numbers. ^'^ The equations of motion 
are integrated with the velocity Verlet algorithm utiliz- 
ing a time step At = 0.005 r. The temperature is set to 
ksT = 1.2 e in all the simulations. At this temperature 
the density of the fiuid is found to be pL = 0.788 
while the vapor density is negligible i^i^i Under these 
conditions, the surface tension is large enough to ensure 
the stability of droplets under an external force. 



III. BULK PROPERTIES 

In order to study the statics and dynamics of droplet 
one needs to evaluate bulk properties of the fiuid as a pre- 
requisite. Quantities, such as surface tension and shear 
viscosity, determine the shape of the droplet and the en- 
ergy dissipation rate, respectively. First, we compute 
the self-diffusion coefficient D of the polymers. The sys- 
tem we use in order to do that is a cubic box of size 
L = 19.04 a containing N = 5440 monomers. The diffu- 
sion coefficient of a polymer is computed from the inte- 
gral of the auto-correlation function of the velocity of its 
center of mass. 



D=- dt (vcM(t) •vcm(O)) 
-J Jo 



(9) 



We find a self-diffusion coefficient D = 0.0157 ± 
0.0030 (T^/r. Moreover, computing the end-to-end dis- 
tance of a polymer we obtain i?ee — \/ (-RL) — 3.447 a. 
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The Rouse relaxation time of the polymer is rp = 
RI^/{^'k'^D) = 25.6 ± 5 tM. This relaxation time per- 
mits to estimate an upper limit of the shear rate that 
one can apply to the fluid before the polymers starts 
to be significantly deformed: The Weissenberg number 
We — jrp should be of the order of unity or less, hence 
7max ~ 0.05 /t. 

The shear viscosity, rj, is obtained from the auto- 
correlation function of the off-diagonal elements of the 
microscopic flux associated to viscosity,— 
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FIG. 2: Sketch of a droplet defining the notations in the text. 



where V is the volume of the fluid and the microscopic 
flux is, 



Ja/3 = XI 



a P 



VPo^p (11) 



1<J 



where P^p denotes the pressure tensor. We use the same 
system parameters as for the calculation of the diffusion 
coefficient to compute the shear viscosity. Averaging over 
2-10^ time units t yields the value ?7 = 5.3±0.1 u"^ j ^pmk. 

Finally, we compute the surface tension, 7, using a slab 
geometry.--'-- In this case the surface tension is found 
from the anisotropy of the pressure tensor according to: 



(12) 



where = = 30 cr are the lengths of the simulation 
box in the y and z directions, respectively. Thus, 2LyLz 
is the surface area of the fluid. The box size is Ly — 
Lz — 30 a and = 100 cr. In order to compute 7 we 
use a larger system than previously, namely = 20 000 
monomers. The value of the surface tension we obtain is 
7 = 0.5154 ± 0.0023 e/cr^. 



IV. STATICS 

In this section, we compute the contact angle and ad- 
hesion energy of droplets for various strengths of the 
polymer-substrate interaction, Cs, and system sizes, A^. 
We use periodic boundaries along the y axis and the 
droplet spans the simulation cell in this direction. As 
a consequence, the droplets are cylindrical. By virtue of 
the cylindrical shape, the radius of the droplet increases 
as ~ \/N (instead of N^^^ for a cap-shaped drop) and, 
thus, larger droplets can be studied. The size of the sim- 
ulation cell along the y axis is fixed to Ly ~ 19.93 cr in all 
the following simulations. One can characterize the ge- 
ometry of a cylindrical droplet by two parameters - the 
equilibrium contact angle, 9e, and the radius, R. We 
have schematically depicted a droplet in Fig. [21 



Other quantities such as volume V, area A of contact 
with the substrate, and height of the center of mass Vz can 
be expressed as a function of R and 9 e ■ For a cylindrical 
droplet one obtains the following relations: 
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(13) 
(14) 

(15) 



We have computed the equilibrium density profiles for 
system sizes ranging from N = 2 000 to A^ = 100 000 
monomers and substrate strengths varying between — 
0.2 e and = 0.8 e. The statistical average of the density 
is typically extended over 20 000 time units r and the 
spatial binning size for the density profiles is chosen as 
1 cr in both directions. The results are depicted in Figs. [3] 
and[31 

The contour line represents the average density, {pl + 
pv)/2, between the density of the dense polymer liquid, 
and the vapor, pv which coexists with it. We note 
that the vapor density is very small (~ 10~^ <^~^)- Uti- 
lizing the relation between center of mass height and con- 
tact angle p5p . one can compute the contact angle in a 
systematic way without referring to the detailed shape 
of the liquid-vapor interface in the vicinity of the three- 
phase contact line. Indeed, one can easily fit a circle to 
the contours of Figs. [3] and S] and obtain the curvature, 
R, of the droplet. On the other hand, the height, r^, of 
the center of mass is straightforwardly monitored in the 
course of the simulation. We depict in Fig. [5] the depen- 
dence of the radius on the number of monomers, A^, and 
the height of the center of mass, r^, as a function of the 
radius for both systems. As we see, linearly increases 
with R showing that there is no significant dependence 
of the contact angle on size. First, the length of the 
three-phase contact line is independent of the droplet size 
for a cylindrical droplet i"^^'^*^ Second, the interface poten- 
tial, which describes the interaction between the liquid- 
vapor interface and the substrate, is rather short-ranged. 
Third, by virtue of the low temperature and the concomi- 
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FIG. 3: Density contours for increasing system sizes varying 
from iV = 2000 to iV = 100000 monomers for ts = 0.4 e. The 
equilibrimn contact angles is 6e = 130°. 

tantly low vapor density and low liquid compressibility, 
the shift of the pressure away from the bulk coexistence 
value due to the Laplace pressure of the curved liquid- 
vapor interface of the droplet does not give rise to a no- 
ticeable change of the liquid or vapor properties. 

The slopes of Figs.[5Ub) and[5jd) yield the equihbrium 
contact angles 9e = 130° and Oe = 29.4° for = 0.8e 
and = 0.4e, respectively. 

We observe that does not vanish for zero radius for 
two reasons: First, the second layer of the substrate is 
approximately located at a height of 0.7 <j. Second, there 
is a depletion layer of the polymer liquid near the sub- 
strate that shifts the height of the center of mass farther 
upwards. We applied this systematic method to other 
substrate strengths in order to evaluate how the contact 
angle and the adhesion free energy depends on eg. In 
equilibrium the macroscopic contact angle of a droplet 
point is determined by the Young equationj^i^ 

isv - ISL - 7 COS 6'b = (16) 

where 751/ , 7sl denote the surface energies between the 
substrate and the vapor phase, and the substrate and 
the liquid phase, respectively. One can rewrite Eq. (jl6p 
in terms of the adhesion free energy per area unit W as, 

W = '^{l + coseE) . (17) 

As a consequence, the fluid will completely wet the 
surface if > 27. When W < 27, however, there will 
be only partial wetting, i.e., droplets with a finite contact 



angle will form. Hence, by representing W as & function 
of the strength of the substrate, £5, one can extrapolate 
towards the value where the wetting transition occurs. 

Alternatively, one can compute the adhesion energy 
from the Molecular Dynamics simulations by computing 
the average energy of the polymer-substrate interactions, 
(Ulj). Utilizing the statistical mechanics definition of 
the free energy 

F^-ksTlnZ where Z ^ J dre-'^^'^% (18) 

denotes the partition function and = e^Vs the 

Lennard-Jones interaction between polymer segments 
and substrate. We diff'erentiate Eq. with respect 
to Eg. This procedure yields dF/deg = (Vb) = (f/£j)/es. 
Integration gives us the free energy as a function of the 
strength of the polymer-substrate interaction: 

F= rde'^^^ + C. (19) 
Jo £ 

where C is an unknown constant of integration. The ad- 
hesion energy W is then simply, W — where A is 
the area of the interface between the fluid and the sub- 
strate. We have computed this adhesion energy for val- 
ues of the strength between eg — 0.2 — 1 e. The fluid is 
confined between two substrates perpendicular to the z 
direction while periodic boundary conditions are applied 
along the x and y directions. The system dimensions are 
— 19.84 a and Ly = 19.93 a while the distance be- 
tween the substrates varies between = 19.74 — 20.96 a 
in order to maintain the coexistence density in the fluid 
sufficiently far away from the substrate (bulk) . The sys- 
tem is comprised of N=5440 monomers. In Fig. [6] we 
present the results obtained by both methods. 

We see that both methods to compute the adhesion 
free energy are in rather good agreement. The constant 
of integration is found to be C « —0.05. This suggest 
the presence of a weak first order or possibly second or- 
der drying transition for small values of Cs in accord with 
previous studies of Lennard-Jones monomer fluids^i^>^ 
and polymeric fluids at a higher temperaturei^i^ The 
small value of the adhesion energy for small values of 
eg allows the droplet to easily detach from the substrate 
due to thermal fluctuations. If an external body force 
is applied, driven droplets will be unstable under these 
conditions. On the other hand, the extrapolation of the 
data for large values of Cs shows that a wetting transition 
of first order occurs at a substrate strength of approxi- 
mately e^* = 0.82 e. 

The fact that is a non-linear function of the sub- 
strate strength, e^, can be rationalized by writing the 
total substrate-fiuid potential energy in the form 

C/tot= / dzp{z)U^^u{z). (20) 
Jo 
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FIG. 4: Density contours for increasing system sizes varying from A'^ = 2000 to A'" = 100000 monomers for es — 0.8 e. The 
equilibrium contact angles is ds = 29.4°. 
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FIG. 5: Radius as a function of the number of monomers and 
the center mass height as a function of the radius for ta = 0.4 e 
in respectively (a) and (b), and for = 0.8 e in (c) and (d). 
The circles are the results of the simulations and the dashed 
line the fits. 



FIG. 6: Adhesion free energy W per area units versus the 
strength of the substrate ta. The circles are the data obtained 
from contact angle computations and with Eq. (|17p while 
the dashed line is from the average Lennard- Jones potential 
energy according to Eq. (|19p . 



While the substrate potential scales like C/waii Cs, the 
density also increases with because the wall becomes 
increasingly attractive giving rise to a strong increase of 
the adhesion energy and the layering at the substrate. 



DYNAMICS 



A. Boundary condition 



dvx 
' dz 



(21) 



that quantifies the balance between the viscous stress due 
to the shear in the fluid and the friction stress at the 
boundary, rj denotes the viscosity of the liquid, and A is 
the friction coefficient, — Vx{zh) denotes the slip ve- 
locity at the hydrodynamic position, of the substrate. 
The quantity S 



When a fluid conflned between two walls is sheared 
with a shear rate 7 a no-slip boundary condition is often 
observed . ^^'^^ However, a finite slippage velocity can oc- 
cur at non-wetting substrateS ]^^i^° The simplest bound- 
ary condition that can describe the fluid velocity at the 
substrate is the Navier slip condition^ 



S = 



A 



(22) 



characterizes the slip length. If the shear rate 7 is not 
too large then the slip length, 6, does not depend on 7 
and it can be computed by non-equilibrium molecular 
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dynamics simulation. To this end, one regards the veloc- 
ity profiles obtained from Couette flow, (i.e., boundary- 
driven planar shear flow) and Poiseuille flow (i.e., a bulk 
force is applied to the fluid). 

The linear velocity profile of Couette flow, 

Vx{z) = ^ [z - Zh - S] (23) 



yields the parameter combination Zh + 5. Additionally, 
using the velocity profile of Poiseuille flow 



This is partially due to the fact that we employ the value 
of the shear viscosity corresponding to the (bulk) coex- 
istence density, pL, of the liquid, while there are impor- 
tant layering effects present at the substrate for large e^, 
enhancing the effective viscosity at the substrate. The 
largest slip length we observe is smaller than 200 a com- 
patible with a recent experiment suggesting that slip 
lengths are of the order of nanometers <^ 

B. Steady state velocity 



Vx[Z) = 

V 



. 1 2 1 2' 

ZhS + -Zf, - -z 



(24) 



one can compute Zh and 5, independently. Alternatively, 
one can also evaluate 5 from equilibrium molecular dy- 
namics by integrating the auto-correlation function of the 
tangential force exerted by the wall on the fluid,-^"^ 



A 



1 



keTA 



dt{F,{t)FM) 



(25) 



and then the slip length is simply given by Eq. where 
•q is the shear viscosity of the fluid in the vicinity of the 
substrate. The position Zh of the interface is about tr, 
which agrees with the shift due to the presence of the 
substrate. In Fig. [7|we show in the slip length computed 
by both the non-equilibrium method and the equilibrium 




FIG. 7: Slip length 5 as a function of the substrate strength 
Es. The plain line is obtained from the Couette and Poiseuille 
profiles while the dashed line from the Green-Kubo relation 
(|25[1 . The inset presents the velocity profiles of the Couette 
and Poiseuille fiow from which the slip length is estimated for 
es = 0.6 e. 



For strong attraction, e^, between polymer and sub- 
strate the Green-Kubo method yields slightly larger val- 
ues of the slip length than the non-equilibrium method. 
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FIG. 8: Sketch of the velocity field inside the droplet. 

In order to study the dynamics of the droplet, we apply 
a constant acceleration a to all particles in the system. 
After a relaxation time, the droplet attains a steady state 
where the energy imparted onto the system by the body 
force is dissipated by the viscous liquid. According to de 
Gennes^ the total energy dissipation rate in the droplet 
can be decomposed in the following way: 



^E=^ E+E+E 



(26) 



TOT 



where is the viscous dissipation due to the rolling 
motion, denotes the dissipation at the three-phase 
contact line, and denotes the dissipation in the pre- 
cursor film. In our system, the third term can be dis- 
carded because there is no precursor film. The second 
term is due to a possible pinning of the particles of the 
fluid to the substrate enhancing dissipation at the con- 
tact line largely. However, since we use a regularly and 
finely corrugated substrate in our study one cam assume 
that such processes do not play an important role. Thus, 
the total dissipation should simply be: 



(27) 



TOT 



However, an extra term should be added to Eq. (P7)) in 
case of a finite slip, ^ > 0, at the substrate. As sketched 
in Fig. [51 the velocity field inside the droplet can be de- 
composed into a rolling motion and the uniform trans- 
lation with the slippage velocity. This slippage at the 



7 



substrate results in a second source of dissipation. Due 
to the corrugation of the substrate, the fluid in its vicin- 
ity undergoes important velocity gradients enhancing the 
dissipation largely. This description is the microscopic 
view of the friction between liquid and substrate. Thus, 
the total dissipation rate can be written as 



(28) 



TOT 



where the term friction dissipation and it is 

proportional to the area of contact A between the droplet 
and the substrate. The second term on the right hand 
side is exactly known: The power dissipated by the sub- 
strate is 



(29) 



where, in the last step, we have approximated the scale 
of the typical shear inside the droplet by U/{S + r^). 

Eq. ((33|) describes the cross-over between two regimes: 
For small drops, SA ^ V or S Vz, the friction at the 
substrate is the dominant dissipation mechanism. In this 
limit, we obtain rjU/pLa ^ V6'^/6A ^ R ^ ^/N. In the 
opposite limit, the viscous dissipation inside the droplet 
dominates and tjU/plci ^ ^ N. This phenomenolog- 
ical argument can be generalized to cap-shaped droplets 
in three dimensions. Then, the rational yields U/ p^a ~ 
R - and U/pLa ^ R'^ ^ iV^/^ for small and large 

droplets, respectively. 

To turn this phenomenological consideration into a 
more quantitative description some information about 
the velocity field inside the droplet is required. The flow 
inside the droplet is governed by the steady-state Navier- 
Stokes equation which in the limit of small Reynolds 
number takes the particularly simple form: 



where the dynamic friction force is 



F, = v^A 
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(30) 



If the droplet is sufficiently flat, the dissipation due to 
the rolling motion can be evaluated in the framework of 
the lubrication approximation. 



tJ2^V / dV 
Jv 



dvx 
dz 



(31) 



where V denotes the volume of the droplet. The viscous 
dissipation scales like the volume V times the square of 
the shear rate. Thus, for flat droplets, one obtains 



(32) 



where ^ is a numerical factor that describes the flow pro- 
file inside the droplet and that only depends on the con- 
tact angle. Here, we assume that the shear at the sub- 
Vs/5 sets the scale of the shear inside 



strate 



dz lz=Zfc 



the entire droplet. This assumption is obviously fulfilled 
for droplets with a small contact angle and a large slip 
length. We will explain below that this viscous volume 
dissipation is most important for flat drops. Moreover, 
we shall verify below that this crude assumption remains 
a reasonably good description even for droplets with a 
larger contact angle (cf. Fig. [9]). 

Now, using the fact that the total energy dissipation 
rate of the system is simply the velocity of the center of 
mass times the bulk force, one arrives at 



2 2 

PLVaU - + ^ 



u 



5 + rz 



{5 A + iV) 



(33) 



(34) 



dz 



PLd 

V 



(35) 



This assumption is justified if the velocity U of the 
droplet is sufficiently small. Thus, the profile of the ve- 
locity Vx parallel to the substrate is a parabolic function 
of the coordinate, z, perpendicular to the substrate^^ 



Vx — az + (iz A- 



(36) 



We additionally assume that the coefficients, a, /3, and 7 
do not depend on the lateral position, x, inside the drop. 
From Eq. [35]we read off the first coefficient, a = — 

Second, we require the tangential stress to vanish at 
the top of the droplet, i.e., at z = H, because the vapor 
has vanishingly small density and viscosity. By virtue of 
the symmetry of the (static) droplet profile, the tangen- 
tial velocity coincides with Vx at the top. 



dvx 
dz 



= 



(37) 



H 



This condition determines the coefficient /3 and one ob- 
tains 



2-1-7 



(38) 



Finally, we utilize the Navier slip condition 
(cf. Eq. (HI])) at the substrate. In the following we 
set the origin of the z axis such that z/i = 0. 



Zh=0 



V 



H (39) 



This equation determines the slip velocity at the sub- 
strate 



Vs = tlb 

V 



(40) 
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and the velocity profile takes the final form 



PLa 



H 



5H 



(41) 



surface friction and viscous volume dissipation regimes 
as 



Eq. (I4ip provides an approximation of the velocity pro- 
file inside the droplet as a function of the droplet shape, 
the slip length, (5, and the force density pLa. 

We emphasize that the height appearing in Eq. ([1T|) 
H = R{1 ~ cos 9e) does not depend on the lateral posi- 
tion, X. This is in marked contrast to previous work(^ 
where H was chosen as the local height h{x) of the 
droplet, which characterizes its shape. This identifi- 
cation, however, would result in a spurious singular- 
ity of the viscous dissipation at the three-phase contact 
linei^i^ Note that the maximal velocity in x direction, 
^max _ Vx\z=h{x): dccrcascs with the lateral distance 
from the center of the drop. This behavior implicitly 
captures some aspects of the rolling motion (cf. Fig.[T|). 

All parameters that dictate the velocity profile in 
Eq. (|4ip have been independently determined: The vis- 
cosity T] and coexistence density pL are known from the 
bulk properties of the fluid, the contact angle 9e from 
static considerations of droplets and the slip length 5 
from the Poiseuille and Couette profiles. One can thus 
predict the velocity profile and quantitatively compare 
them to the results of the molecular dynamics without 
any adjustable parameter. This comparison is presented 
in Fig. El 
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FIG. 9: Velocity profile as a function of the height z for dif- 
ferent lateral positions, x, in the droplet. The inset shows the 
sections which have been used to calculate the velocity profile 

V:,:{X, z). 

The comparison validates the basic assumption that 
the lateral velocity Vx is independent of the lateral po- 
sition, X. Indeed, the simple analytical expression (j4ip 
provides an excellent description of the simulation data 
using the independently determined parameters. 

Considering the velocity at the top of the droplet — 
Vx{H) one can evaluate the crossover radius Rc between 



Rc 



25 



1 — cosOe 



(42) 



Thus in order to have a large crossover radius, Rc, one 
should have a large slip length 5 and/or a small contact 
angle, Oe- 

The detailed velocity profile inside the droplet may be 
difficult to measure experimentally but the cross-over be- 
tween the two regimes can be observed by the steady- 
state velocity of the center of mass. Approximating the 
center of mass velocity simply by the velocity Vx{rz) at 
the height of the center of mass r^, we obtain: 



U = 



PLa 
V 



SH 



(43) 



In order to test our model we computed the depen- 
dence of the steady-state velocity on the droplet size for 
two substrate strengths, = 0.4 e and eg = 0.8 e. We 
first note that our model is valid only for small Reynolds 
numbers such that no inertial effects are present and the 
droplet shape resembles its equilibrium shape. Moreover, 
in this linear regime the center of mass velocity is pro- 
portional to the acceleration and, thus, we plot the ratio 
U/a instead of U. 

Care must be taken in choosing the acceleration, a. 
A too small value of a gives rise to impractically long 
computation times and large statistical errors in the ve- 
locity profiles and the average center of mass velocity. 
If one chooses a too large, however, non-linear effects 
become important and one tends to underestimate the 
ratio U /a. Since at a fixed acceleration, a, the velocity 
of center of mass increases with size, one reaches the non- 
linear regime when the droplet size grows. We note that 
at very high velocities two different behavior can be ob- 
served: For small contact angles, the droplets break up 
into smaller pieces, while for large contact angles they 
detach from the substrate. 

In Fig. [TO] and Tab. U we present the results of the sim- 
ulations and compared them to the predicted behavior 
. We observe that Eq. (|43p is in excellent agreement 
with the results of the Molecular Dynamic simulations. 
The error bar on the prediction characterize the uncer- 
tainty associated with the error of the input parameters. 
The figure also depicts the asymptotic limits for small 
and large droplets. Our simulation data are clearly in 
the cross-over region. The data for the droplet with the 
larger contact angle (cf. panel a) are mainly in the sur- 
face friction dominated regime, while the flatter droplets 
(shown in panel b) already reach the regime of volume 
dissipation. The crossover between surface dominated 
dissipation and volume dominated dissipation occurs at 
Nc ~ 160 000 and « 33 000 monomers for = 0.4 e 
and Eg = 0.8 e, respectively. Consequently, surface ef- 
fects are more important than volume dissipation in the 
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FIG. 10: Velocity ol the center of mass per unit acceleration 
as a function of the number of monomers for (a) = 0.4 e 
and (b) es = 0.8 e. The circles are the results of the molecular 
dynamics simulations and the dashed lines the analytical re- 
sults. The error bars mark the uncertainties associated with 
the independently determined input parameters. The asymp- 
totic limits for small and large droplets are indicated by the 
solid lines proportional to y/N and N , respectively. 





iV 


Ur/a 


ar^ /a 





4 


2000 


0.0056 


0.0001 





4 


5000 


0.0114 


0.0001 





4 


10000 


0.0180 


0.0001 





4 


20000 


0.0272 


0.0001 





4 


50000 


0.0230 


0.00005 





4 


100000 


0.0364 


0.00005 





4 


200000 


0.0295 


0.000025 





8 


2000 


0.000367 


0.0001 





8 


5000 


0.000699 


0.0001 





8 


10000 


0.00114 


0.0001 





8 


20000 


0.00185 


0.0001 





8 


50000 


0.00185 


0.00005 





8 


100000 


0.00162 


0.000025 





8 


200000 


0.0027 


0.000025 



TABLE I: Parameters and results of the Molecular Dynamics 
simulations of droplets driven by an external body force. 



droplet with the large contact angle. This can be ratio- 
naUzed as follows: In droplets with large contact angles 
the velocity profile in the upper part of the droplets is 
close to a purely rotational motion which is characterized 
by very low dissipation.'^ This is clearly seen in Fig. [1] 
which depicts the velocity profile in the comoving frame 
for a droplet with N = 100 000 monomers and a bulk 
acceleration a — 0.00005 cr/r^. Thus, the viscous dissi- 
pation in the inside of the droplet is small although the 
volume is large. 

This explanation is corroborated by the profile of the 
viscous energy dissipated per unit volume of the fluid. In 
the inset of Fig. [T] we show this quantity as a function 



of the height, z. Unfortunately, we are unable to calcu- 
late the analog quantity for the droplet with the smaller 
contact angle with our computational resources because 
for flat droplets the linear regime is limited to very small 
accelerations and the velocity of the center of mass is 
orders of magnitudes smaller than the thermal velocity, 
VT = ^kBTjm = 1.095 cr/r. 

As we see in the Fig. [T] the energy dissipation rate 
in a droplet with a large contact angle is twice as large 
in the lower region than it is in the upper one where 
the stream lines resemble those of a purely rotational 
flow. On the other hand, if we decrease the contact angle, 
the flow inside the droplet will increasingly deviate from 
this Huygens motion and, consequently the region where 
important viscous dissipation occurs will expand. As a 
result, for a small droplet of fixed volume, if the area 
of contact with the substrate is large (Ge small) volume 
dissipation dominates, while if the area of contact is small 
(Qe large) surface dissipation dominates. 



VI. CONCLUSION 

In this paper we studied equilibrium and non- 
equilibrium properties of polymer droplets in contact 
with a corrugated substrate by Molecular Dynamics sim- 
ulations and quantitatively compared the results to a 
simple analytical model. First, we measured the contact 
angles for different substrate strength via a geometrical 
method that does not rely on the detailed shape of the 
liquid-vapor profile at the three-phase contact line. Ad- 
ditionally, we evaluated the adhesion free energy of the 
droplet and found good agreement between the shape of 
the cylindrical droplets and the macroscopic contact an- 
gle predicted by the Young equation for our coarse-grain 
polymer model. In addition to these static characteristics 
we studied non-equilibrium properties of the system such 
as the boundary condition (i.e., hydrodynamic interface 
position and slip length) and the steady-state velocity of 
the droplet when a bulk force is exerted on it. 

An analytical relation predicting the velocity field in 
the direction of the flow is derived and the scaling of the 
velocity of the droplet's center of mass discussed. We 
show that there are two different regime in the droplet 
dynamics: For small droplets, friction at the substrate 
dominates dissipation and the droplets velocity, J7, in- 
creases like \/]V with the number of particles, N ^ at con- 
stant acceleration. For large droplets, however, the vis- 
cous dissipation of the flow in the interior of the droplet 
dictates the behavior and one observes [/ ~ iV. Surpris- 
ingly, we find that the larger is the contact angle, and 
thus the smaller is the area of contact, the larger are the 
surface effects. 

While our study is restricted to cylindrical droplets 
for computational constraints the qualitative arguments 
carry over to cap-shaped droplets. The cross-over be- 
tween the surface- and the volume-dominated dissipa- 
tion behavior occurs at a droplet radius on the order 
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Rc — cos9e)- Using a slip length on the order 

of micrometer as extracted from recent dewetting ex- 
periments by Fetzer at al~^ on low-molecular-weight 
polystyrene melts on octadecyl- (OTS) and dodecyl- 
trichlorosilane (DTS) polymer brushes and a contact an- 
gle of Qe ~ 67° we predict Rc ~ 3/im which is in an 
experimentally accessible range. 
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